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ABSTRACT 

A new model for the high-energy emission from pulsars is developed by considering 
charged particle motion in the fields of a spinning, highly magnetised and conducting 
sphere in vacuum. A generally applicable approximation to the particle motion in 
strong fields is developed and applied to the numerical modelling, and the radiation 
emitted by curvature emission is summed to generate light curves. The model predicts 
many of the observed features of pulsar light curves. This paper outlines the basic 
properties of the model; a subsequent paper will discuss the statistical properties of a 
population of model pulsars and apply the model to the known gamma-ray pulsars. 
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1 INTRODUCTION 

Much work has been done over the past twenty-five years 
on the low- frequency radio emission detected from pulsars, 
in large part due to the huge volume of data available. In 
the last decade, however, attention has been turning toward 
the high end of the emission spectrum, as more reliable data 
become available (see, for example, Grenier et al 1993 and 
Masnou et al 1994). 

The main problem in developing a theoretical model 
of the pulsar magnetosphere is the complexity of finding 
a self-consistent solution to the set of equations defining 
the electromagnetic fields and particle densities and veloci- 
ties as functions of time. Theoretical attemps to model the 
high-energy emission generally make assumptions about the 
global structure of the magnetosphere, and then restrict 
their attention to local sites where large electric fields par- 
allel to the magnetic field lines accelerate charges. The two 
most developed examples are the polar cap model of Daugh- 
erty & Harding (1994) and the outer gap model created by 
Cheng, Ho, & Ruderman (1986) and further developed by 
Chiang & Romani (1994). 

We make a set of assumptions about the global magne- 
tospheric structure that differs from past models: we assume 
that the charge density is very small in the outer magneto- 
sphere, and that the pulsar can be treated as a spinning, 
highly magnetised and conducting sphere in vacuum (the 
electromagnetic fields around such a star were derived by 
Deutsch 1955). This assumption is justified in our model 
because we find that any charges created in the inner mag- 
netosphere follow paths constrained to stay close to the star. 
They are thus unable to populate the outer magnetosphere 



and short out the electric fields (the distinction betwee n th e 
inner and outer magnetosphere is elaborated in section 3.2). 



By so fixing the fields, particle motion can be integrated 
relatively easily, and the radiation observed from the model 
pulsar can be estimated fairly accurately. With low number 
densities, inverse Compton scattering and photon/photon 
interactions are unlikely, and synchrotron emission is not 
important, so curvature radiation is the only source of emis- 
sion. The back reaction on the particles of the curvature 
emission is included in the determining the particle ener- 
gies. 

An approximation to the particle motion in strong fields 
in developed which is applicable in all realistic pulsar mag- 
netospheres where gamma-ray production is efficient. The 
approximation is used to simplify our numerical integration, 
but is general enough that it may find application in other 
models of the pulsar magnetosphere. 

In this paper (Paper I), the basic characteristics of the 
model will be examined. Higgins & Henriksen (1996) (here- 
after Paper II) will apply the model to the known high- 
energy pulsars, consider how the model pulsar ages, and 
estimate the number of high-energy pulsars that should be 
found in gamma-ray surveys. 



2 THE DEUTSCH FIELDS AND THE 
VACUUM APPROXIMATION 

2.1 The Deutsch Fields 

By fixing the electromagnetic fields, the problem of calcu- 
lating the radiation pattern changes from a system of par- 



© 1996 RAS 



2 M.G. Higgins & R.N. Henriksen 




Figure 1. Selected closed field lines for the orthogonal Deutsch 
field pulsar in the plane perpendicular to the rotation axis. The 
circle shows the light cylinder, and the dimensions are in units of 
the light cylinder radius R\ c . 



tial differential equations to a system of ordinary differential 
equations, and becomes much more tractable. 

The Deutsch fields look like the fields of a spinning mag- 
netic point dipole at radii much larger than the stellar ra- 
dius, with wrapped-up magnetic field lines outside the light 
cylinder and large components of electric field parallel to the 
magnetic field. Inside the star the electric field is frozen-in, 
and this boundary condition at the stellar surface signifi- 
cantly perturbs the fields near the star. Field lines in the 
plane perpendicular to the rotation axis for the orthogonal 
Deutsch fields (when the magnetic dipole moment is per- 
pendicular to the rotation axis) near the light cylinder are 
shown in Figure [l] for selected closed field lines and Figure ^ 
for open field lines. 

Four parameters entirely characterise these fields: the 
angular velocity of the star fl, the radius of the star a, the 
magnitude of the effective magnetic dipole moment n, and 
the inclination \ of the dipole moment with respect to the 
rotation axis. 

An estimate of the electric and magnetic field strength 
is useful for the calculation of various quantities; the order of 
magnitude of the electric and magnetic fields near the light 
cylinder is approximately 



E~B 



r3 



(1) 



The vacuum approximation holds as long as the charge 
densities and current densities give negligible perturbations 
to the fields. This is true when the charge density is much 
less than the Goldreich-Julian charge density: 
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Figure 2. Selected open field lines for the orthogonal Deutsch 
field pulsar in the plane perpendicular to the rotation axis. The 
circle shows the light cylinder, and the dimensions are in units of 
the light cylinder radius R\ c . 



where the fields are frozen-in. Here, Q is the angular ve- 
locity of the star and B is the magnetic field at the point 
in question. This corresponds approximately to a number 
density: 



tigj — 4 x 10 4 cm 



;) 4 ( 



100rads/s' v 10 30 G-cm 



(3) 



where Ri c = c/fi is the radius of the light cylinder. Note 
the strong variation with the angular velocity of the star; a 
relatively small charge density at the light cylinder of a slow 
pulsar can alter the fields. 

The fields far from the star look like a spherically- 
propagating electromagnetic wave with frequency equal to 
the rotational frequency of the star. This radiation removes 
rotational energy from the star, slowing it with time. There- 
fore, the effective magnetic dipole moment of the Deutsch 
fields can be estimated from the period and period derivative 
of the pulsar rotation: 



fi ~ 3.2 x 10 37 (PP) 



( i 

37/ n ml/2 1 10 45 g-cm 2 
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(4) 



The moment of inertia of neutron stars is unknown by 
up to a factor of 10 depending on the equation of state; 
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this introduces some uncertainty into the calculation of the 
effective magnetic dipole moment. 

This result is similar to that of magnetodipole radiation 
from a point dipole derived by Ostriker & Gunn (1969), since 
the Deutsch fields only differ from the point dipole fields near 
the star. 



2.2 Particle-Photon and Photon-Photon 
Interactions 

This model assumes that the only emission contributing to 
the light curves is a result of curvature radiation by the high- 
energy particles moving along curved paths in the magneto- 
sphere. This is justified as long as the number densities are 
small; how small can be roughly calculated by examining the 
mean free path of photons subject to Compton scattering off 
electrons and to self-scattering. 



2.2.1 Inverse Compton Scattering 

The cross-section for Compton scattering (or inverse Comp- 
ton scattering) is approximately equal to the Thompson 
scattering cross-section for photon energies in the electron 
rest frame much less than the electron rest mass. For photon 
energies in this frame much greater than the electron rest 
mass, the cross-section is approximately jarmc 2 / E-y (Jack- 
son (1975)). Therefore, the Thompson cross-section will be 
used as an upper limit to the cross-section (although this ne- 
glects resonant interactions with the magnetic field, which 
can increase the cross-section by several orders of magnitude 
above the Thompson value for a limited energy range - see, 
for example, Daugherty & Harding 1991). 

The optical depth of the magnetosphere to inverse 
Compton scattering for a photon is then approximately 



TICS 



: n c <TT-Ri< 



where n c is an average electron density, <tt = 6.65 x 10 -25 
cm 2 is the Thompson scattering cross-section, and R\ c = 
c/f2 is the light cylinder radius. 

The number density of photons in the magnetosphere 
can be roughly estimated from the power output of each 
electron P e , the average photon energy i? 7 , and the charac- 
teristic time the photons are within the magnetosphere 1/Q 
as 



(5) 



SI E-y 



To ensure that the number of photons which scatter off 
electrons is much less than the number of electrons, 

n 7 ncs < «. e 

This puts a limit on the number density of electrons of 

0,^ E-y 



CTtcPo 



(6) 



For relativistic motion on a path with radius of curva- 
ture p, the power and average frequency of curvature emis- 
sion can be written (e.g., Jackson 1975) 



2e 2 c 4 

"e = 0-5-7 
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p 



where 7 is the charge's Lorentz factor and e its charge. 
Substituting these into inequality (H) , 
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(9) 



The Lorentz factor of particles can be estimated from 
the fields (see equation (|l4[)) by equating the power output 
from radiation with the power input from the electric field. 
Roughly, this gives 



7 



This can be used to define an upper limit to a com- 
bination of the stellar parameters below which this model 
remains applicable by setting the Golreich- Julian number 
density equal to the limit of equation (|^) (the number den- 
sities must always be much less than the Golreich- Juliann 
number density). Using the expression for 7 above, the con- 
dition for neglible scattering is 



fi^M 3 < 3.3 x 10 4 



(10) 



For example, with a magnetic dipole moment of 10 G 
cm 3 (giving a surface magnetic field of approximately 10 12 
G), the upper limit on SI is approximately 5000 rad/s, much 
larger than any standard pulsar spin rates. Thus, for the 
population of normal pulsars, the Goldreich- Julian number 
density is not sufficient to ensure significant Compton scat- 
tering. 

2.2.2 Photon-Photon Scattering 

A similar calculation can be done for the case of photon- 
photon scattering, where the cross-section above the pair 
production limit is again maximally of order the Thompson 
cross-section. For the spectrum to remain unperturbed as it 
leaves the magnetosphere, 

t 7 _ 7 <C 1 

where r 7 _ 7 = fiphotonCT-Ric- This gives exactly the same 
constraint on the electron number density and on the angu- 
lar velocity and dipole moment (equation die])) as the inverse 
Compton scattering criterion above. 

While the fraction of photons affected by photon- 
photon scattering is rather small, the photons which do 
interact create charges in the magnetosphere. This could 
provide a source for the charges, though a more complex 
calculation is required to properly examine the viability of 
this option. 



2.3 Synchrotron Radiation 



A key approximation of the model (see section 3.1) is that 



charges move parallel to the magnetic field line in a frame 
where the electric and magnetic fields are parallel. There- 
fore, in the framework of the model, there is no contribution 
at all from synchrotron radiation. The results of simulating 
particle motion under this approximation have been com- 
pared to those found in a full-blown simulation of the par- 
ticle motion where the approximation is relaxed. The two 
calculations agree remarkably well whenever the fields and 
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c f3E z „ cB — (3E y 
7 [7 j U 



(13) 



Figure 3. A boost along the x-axis shown in the figure moves into 
a frame where the electric field E' is parallel to the magnetic field 
B' . In the pulsar rest frame, the z-axis is taken to be parallel to 
B, and E lies in the y-z plane. In the boosted frame, the charges 
are assumed to move either parallel or anti-parallel to B' . 



energies are large enough to produce significant amounts of 
radiation. 



3 PARTICLE MOTION AND EMISSION 

3.1 The Drift Frame Bead-on-a-Wire 
Approximation of Particle Motion 

The numerical integration of the charged particle motion in 
the simulations is simplified enormously under the applica- 
tion of a key assumption in this work, labelled the drift frame 
bead-on-a-wire (DFB) approximation. We assume that, in a 
frame where the electric field is parallel to the magnetic field, 
the charge has a velocity directed parallel or antiparallel to 
B with magnitude approximately equal to c; the velocity in 
the pulsar rest frame at any point can then be determined 
by boosting back to the original frame of reference. 

Consider a point in the magnetosphere with electric 
field E and magnetic field B. Define the z-axis to be parallel 
to B, and E to have components only in the z and y direc- 
tions. There is a boost of v = cf3x which defines a frame of 
reference where E' is parallel to B' as long as the electric 
field is not perpendicular to the magnetic field and equal 
in magnitude (primed quantities are those in the boosted 
frame). This is shown in Figure |§| the boost can be shown 
to be 



B 2 + E 2 — y/(B 2 + E 2 ) 2 - 4B 2 E 2 



IB Ey 



(11) 



where all the quantities relate to the pulsar rest frame. This 
is valid for E > B or E < B, and even E — B as long as 
the two are not perpendicular. 

We assume that, in this frame, the velocity of the charge 
is either parallel or anti-parallel to B' (depending on the 
sign of qE ■ B, where q is the particle's charge) and that the 
particle is ultra-relativistic, so its speed is nearly c. Boosting 
back to the initial frame, the charge's velocity is then 



cB' 



cpx ± — — -y ± — —y 

7 B' j B' 



(12) 



where U = ^ B 2 + /3 2 E 2 - 2f3E y B and 7 = l/^/l-fi 2 . 
The expression is exact, given the assumption about the 
particle motion in the boosted frame of reference. 

Given an initial position, the charge's motion is now 
defined by the solution to a set of three first-order ordinary 
differential equations if = df/dt, which can be solved numer- 
ically with relative ease. The particle's energy can be found 
by integrating the power input due to the electric field and 
the power output due to the radiation loss. The simulations 
of particle motion under this approximation have been com- 
pared to those when the DFB approximation is relaxed in 
a more sophisticated (and time-consuming) simulation and 
are found to match very closely for all the known gamma- 
ray pulsars (the constraint on the applicability of the DFB 
approximation is considered below). 

It is generally found that, since the fields are so large, 
the radiation reaction force acts like a thermostat, balanc- 
ing the energy input from the electric force with the energy 
output due to radiation. This gives a convenient description 
of the particle's energy at any point, 

where the radiative power output was taken as equation (0). 

The whole approximation is only valid where the 
timescale to radiate away the momentum perpendicular to 
the field direction is much less than the timescale 1/S1 over 
which the fields vary significantly; this reduces to another 
constraint on Q and fj,: 



0,—^ > 5 x 10 7 



(15) 



For a typical fx = 10 30 G cm 3 , this gives S7 3> 6 rad/s, 
which is true for all the observed gamma-ray pulsars. 



3.2 Types of Particle Motion and Justification of 
the Vacuum Assumption 

Three different types of particle motion were observed in 
the simulations. Some charges follow open paths and escape 
to infinity, radiating quite powerfully inside the magneto- 
sphere; these are labelled "outgoing" paths. The second and 
third types of path are closed within the inner magneto- 
sphere. Charges on the second type of path follow a course 
which directs them into the star; on the third variety, the 
charge paths attempt to cross a surface where E ■ B goes to 
zero. Particles encountering such a "null surface" lose their 
energy almost immediately since the radiation power output 
is no longer balanced by an energy input from the electric 
field. Any charge encountering a null surface is forced to re- 
main on it; the electric field switches sign on either side of 
the surface, effectively trapping the particle. They can travel 
out of the magnetosphere on these surfaces at relatively 
low energies without contributing to the observed radiation, 
but perhaps establishing the required charge balance for the 
steady-state magnetosphere. These surfaces then are not so 
much "gaps" as "conduits" for low-energy plasma. 

Which type of path a charge will follow depends on its 
starting position. The fraction of paths which are outgoing 
drops to zero at a radius of approximately 0.1i?i c ; all the 
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charges inside this radius are bound to the inner magneto- 
sphere, except those on the surface with E ■ B = 0, which 
can escape the magnetosphere but at very low energies. This 
defines the boundary between the inner and outer magne- 
tosphere: inside the inner magnetosphere, all particle paths 
are bound. 

This constitutes our justification of the evacuated na- 
ture of the magnetosphere; charges created (presumably due 
to 7 — B pair creation near the star) will not be able to pop- 
ulate the outer magnetosphere, since they are locked into 
bound orbits. We take the bound charges into account by 
increasing the effective stellar radius to 0.1i?i c . In doing so, 
the charges remain bound, even though the field structure 
changes to become "frozen-in". 

When charges are started relatively far inside the mag- 
netosphere, but outside 0. li?i c , most of the paths are still 
bound, and the outgoing paths extend from a small region 
near one of the poles (which one depends on the sign of the 
charge of the particle). 



4 OUTLINE OF THE MODEL 

With this justification of the vacuum assumption for the 
magnetosphere, the basic model can be defined as a rough 
approximation to a global self-consistent solution to the pul- 
sar problem. The outer magnetosphere (r > 0.1i?i c ) is the 
region of particle acceleration, where charges can be ejected 
from the magnetosphere. The only radiation which is impor- 
tant is curvature radiation, for the reasons discussed earlier. 

Particles radiate gamma-ray photons, some toward the 
star, where they create electron-positron pairs on the pow- 
erful magnetic fields in the inner magnetosphere. These 
charges are all bound on closed paths to the inner mag- 
netosphere, and cannot populate the outer magnetosphere 
to short out the large electric fields. A charge balance is es- 
tablished by having charges being created near the star and 
flowing out with low energies along the null surfaces. The 
effect of the outflowing charges on the fields is neglected. 

Since the paths of the charges are all bound for radii 
r < O.f-Ric, charges were started on a sphere of radius 0.2i?i c 
(called the "starting sphere") at intervals of n/16 radians in 
the toroidal and poloidal angles to fully sample the range of 
initial starting positions. The starting energies were unim- 
portant, since the charges immediately move to balance ra- 
diation power out with electric field power in. The source 
for these charges is left undetermined. 

The radiation from these trajectories was calculated and 
binned to build up light curves and spectra for different 
viewing directions; this is discussed in detail in sections ^ 
and ^ 

To simplify the numerical calculation, we took advan- 
tage of three symmetries of the problem. First, there is a 
symmetry (Symmetry 1) in the fields E — * —E and B — > — B 
under a translation — > <^> + tt, 6 ^ n — 8, where 9 is the 
poloidal angle and <f> is the toroidal angle relative to the 
rotation axis. Therefore, if a charge q started at fo follows 
the path f(t), the same charge started at — fo will take the 
path — f(t). The symmetry applies for any orientation of the 
dipole moment to the rotation axis. 

The second symmetry (Symmetry 2) is between charges 
of different sign, but only applies when the dipole axis is 



perpendicular to the rotation axis. The path of a charge — q 
started at (ro, 6, 0) is equal to the path of a charge q started 
at (ro, 8, (j> + 7r) rotated through 180 degrees toroidally. 

The third symmetry (Symmetry 3) is a result of another 
property of the fields, that E(r,d, cj>,t) equals E{r,6,<j> — 
fit, 0) rotated through an angle fit in <j>. This means that 
a charge q started at (ro,9,cf>,to) will follow the path of a 
charge q started at (ro,8,(j> — fito, 0) rotated through fito 
toroidally. This symmetry holds for any orientation of the 
dipole moment. 

Symmetry 1 means than runs need be done only for 
particle starting points in the top half of the magnetosphere; 
paths in the bottom half are reflections of those in the top 
half. A consequence of Symmetry 3 is that charges started 
at a time to not equal to zero can be found from the paths 
started at to — by a rotation azimuthally. Symmetry 2 
means that positron paths can be inferred from electron 
paths, but only in the orthogonal Deutsch fields. 

The combination of Symmetry 2 and Symmetry 3 has 
an important consequence for the inferred light curves when 
the dipole moment is perpendicular to the rotation axis: 
positrons will radiate in exactly the same manner as elec- 
trons, but a half period later. If the electron contribution is 
a single peak, the positrons will contribute a second peak, 
exactly 0.5 later in phase, if the distribution of positrons is 
equal to that of electrons in the magnetosphere. This does 
not hold for non-orthogonal systems, as Symmetry 2 is bro- 
ken. 



5 LIGHT CURVE GENERATION 

With paths of all the charges known, a pattern of radiative 
emission for different retarded times r (r = t — n b s • r/c, 
with n bs a unit vector pointing toward the observer) could 
be generated. 

A numerical routine stepped along the paths for each 
of the different starting angles, and for 100 steps in starting 
angle to- The emission from paths in the bottom half of 
the magnetosphere (and for positron paths when considering 
orthogonal fields) was included by inferring their paths from 
the three symmetries discussed earlier. 

Since the emission generated is only for a single elec- 
tron started at (9, 4>) on the starting sphere and at time 
to, the power must be scaled by a distribution function / e , p 
of charges over the starting sphere and over starting time, 
where f C:P dS}dto is the number of charges (either electrons 
or positrons) flowing through the starting sphere between 
solid angle fi and fi + dfi at times between to and to + dto- 
This distribution function can be varied to model different 
sources; however, since most emission comes from the out- 
going paths, which only originate on a small region of the 
starting sphere, changes in / e , p do not much affect the light 
curves. Throughout this work, f e = f p = f for the light 
curve generation; the actual value of / is constrained by 
the total power emitted by the pulsar, and sets the number 
density at the starting sphere. 

The spectrum of the radiation at one point along a path, 
for a given particle Lorentz factor and path radius of curva- 
ture, is assumed to be instanteously that for circular motion; 
this is discussed in section \l\ The emitted radiation was cal- 
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Figure 4. Representative light curves produced by the model for different viewing directions, spin rates, field strengths, and dipole 
inclinations. See the text for a complete discussion. Time is plotted along the horizontal axis, and arbitrary intensity along the vertical. 
All plots have radio peaks at r = and/or r = 7r/f2, presuming the radio emission to be generated close to the star and near the pole. 
Plot (d) is shown offset in time to show the single peak more clearly. 



culated at several different frequencies to capture the spec- 
trum at each time on the light curve. 

The emitted power was binned into angular bins of 7r/16 
radians in azimuthal and polodial viewing angles, bins of 
T/20 in t (where T is the pulsar period), and bins of 0.25 in 
logo;, where hu> is the photon energy. Any emission which 
was beamed through the inner magnetosphere (r < 0.1i?i c ) 
was not included in the light curves, as it was assumed to 
have been absorbed in the denser plasma near the star. 

Figure ^] shows some representative 1 MeV light curves 
generated by the model, which cover many of the variations 
seen in the observed gamma-ray pulsar light curves. Plots 
(a) and (b) correspond to a Geminga-like pulsar, with an 
angular velocity of SI = 26.5 rad/s and a magnetic dipole 
moment of fi — 10 30 G cm 3 aligned perpendicular to the 
rotation axis. The first is for a viewing angle of 9 V = 84.4°, 
and the second represents radiation into V = 50.6°. These 
show the characteric symmetry between the first and second 
halves of the light curve due to the action of Symmetry 2 in 



the orthogonal field. Light curve (c) is for the same Geminga- 
like pulsar, but with the dipole moment inclined at 45° to 
the rotation axis, seen at a viewing angle of 73°. The two 
main peaks are separated by a phase of 0.45. (d), (e), and 
(f) all correspond to emission from a Vela-like pulsar, with 
fi = 70.6 rad/s and fj, = 3.7 X 10 30 G cm 3 . The dipole 
moment in (d) and (e) is nearly aligned at only 11.25° to 
the rotation axis; (f) corresponds to orthogonal fields, (d) 
represents emission into a viewing angle of 16.9°, and shows 
a well-defined single peak, (e) is seen from a viewing angle 
of 61.9° and shows three quite well-defined peaks (there is a 
suggestion of a triply-peaked light curve for PSR B1706-44 
- Thompson et al, 1995). (f) is seen from a viewing angle of 
39.4°, and shows no discernable pulsation at frequencies less 
than 50 GeV. This pulsar would likely be missed in a pulsar 
survey. 

All the plots discussed relate to emission beamed into 
4> v — 0; one of the poles points toward the observer at t = 
at this angle. The light curves for 4> v 7^ are identical except 
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for a phase change. For all these plots, the radio emission 
peaks (presumed to come from the poles near the star) would 
be seen at r ~ and/or r = n/£l; the gamma-ray peaks are 
significantly displaced in some cases, while in others they 
align quite closely. 



6 TOTAL EMITTED POWER 

The total power emitted into all directions by the electrons 
or positrons at any time t is given by: 



fe, P (*)= / / [f e P e + f P P P ]dQdt 



(16) 



where P e . 



, to; t — to) is the power output of one electron 



(e) or positron (p) at time t — to along the path, started at 
(9,4>) on the starting sphere at time to. / e , p are the distri- 
bution functions of charges on the starting sphere, which is 
not fixed in the model. 



6.1 The Required Number Density 

If the total high-energy power output of a pulsar is known, 
this fact can be used to estimate the number density of 
charges at the starting sphere required to supply that en- 
ergy. If this number density is larger than the Goldreich- 
Julian number density, the fields will be altered and the 
model will break down. 

Taking a constant distribution f c = f p = f, the total 
energy radiated in one period is given by 

E = f / /(P e + P p )dtdt dtt (17) 
J n J t J t 

Since the radiation reaction force fixes the radiated 
power to be very closely equal to the power input from the 
electric field, 



E ■ 0e,p 



(18) 



where £ is a dimensionless electric field which no longer 
depends on the values of /x and f2, and /3 e , P = Se,p/c. 

Defining to = fito and t = Qt, the total energy emitted 
in a period can now be written 



■Ctot — 7 — 5~ £ 



(19) 



where £ is a dimensionless function of the dipole inclination 
angle \ alone: 

£= / E ■ e + P P )dQdT dT (20) 

For example, £ = 36 for the orthogonal fields, and 6.8 
for x = 11-25°. 

Since the total gamma-ray power output from the model 
must be constant with time, the total energy emitted in a 
period is equal to P Q b s 4jf , where P t, s is the observed gamma- 
ray power. 

The value of / can then be related to the observed power 
output: 



2tt c 2 P obs 



eQ. 2 fi£ 
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Figure 5. Distribution of emitted power at 1 MeV with viewing 
angle for two Vela-like pulsars. The solid line shows the power 
emitted by a pulsar with orthogonal fields, and the dashed line 
shows the corresponding curve for a pulsar with a dipole axis 
inclined at 11.25° to the rotation axis. The axes show viewing 
angle in degrees on the horizontal and emitted power scaled to 
unity at the maxima on the vertical. 



The value of / sets the number density at the starting 
sphere required to produce the observed power output, since 
/ = 47rc(0.27?ic)V Therefore, 

The number density can be described somewhat more 
meaningfully as a fraction of the Goldreich-Julian number 
density (equation (^)) at the starting sphere (£ is taken to 
be equal to 36, the value for the orthogonal fields). 



n 
riGj 



' 4.7x10" 



5 cm" 3 ( 



:)( 



r 4 (- 



10 34 ergs/s M 100rad/s' v 10 30 Gcm 3 

This quantity is tabulated in Table 1 for the known 
gamma-ray pulsars, and is less than unity, and therefore 
consistent with the model, in all cases. The values of /x were 
calculated with equation (ft) taking sin x = 1 (and are there- 
fore a lower limit). 

The most serious problem with this model is its inability 
to predict the total power output; instead, it is only shown 
that the observed power outputs require number densities 
less than the Golreich- Julian number density, and therefore 
are consistent with the assumption of vacuum around the 
star. The source of the charge on the starting sphere is left 
unexplained. 



6.2 The Gamma-Ray Efficiency 

The gamma-ray efficiency of a pulsar is defined as the ratio 
of the power output in gamma-rays to the total energy loss. 
For the Deutsch field pulsar, the rate of rotational energy 
loss is 



2n 4 /i 2 sin 2 x 
2.5 x 10~ 32 fi 4 /z 2 



(24) 
(25) 



and the gamma-ray power emitted in our model is given by 
2£ 

P-y = 25 ec / in ( 26 ) 
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Pulsar 


P 


P 


ft 


-Pot 


ra/n GJ 




Reference 


GGmingcL 


0.237 


11.0 


1.63 


1.6 x 10 34 


5.6 x 10" 


3 


Mayer-Hasselwander et al 1994 


Vela 


0.0893 


125 


3.4 


6.3 x 10 33 


1.0 x 10' 


5 


Kanbach et al 1994 


Crab 


0.0333 


421 


3.8 


10 36 


2.6 x 10" 


5 


Ulmer et al 1995 


PSR1509-58 


0.150 


1540 


15 


6.3 x 10 35 


4.2 x 10" 


4 


Laurent et al 1994 


PSR1706-44 


0.102 


93.0 


3.1 


6.1 x 10 34 


2.1 x 10" 


4 


Thompson et al 1992 


PSR1055-52 


0.197 


5.8 


1.1 


9.4 x 10 33 


3.5 x 10" 


3 


Ficrro et al 1993 


PSRB1951+32 


0.0395 


5.85 


0.49 


1.3 x 10 3 4 


4.0 x 10" 


5 


Ramamanamurthy et al 1995 



Table 1 . Number densities at the starting sphere required to give the observed power outputs for the known 
gamma-ray pulsars, as a fraction of the Goldreich-Julian number density. Periods P are given in seconds, 
period derivatives P in 10~ 15 s/s, dipole moments fi in units of 10 30 G cm 3 , and gamma-ray power outputs 
Ptot in ergs/s. All known high-energy pulsars have n/nQj <C 1- 




Figure 6. The variation of f2 p , the fractional solid angle into 
which radiation is beamed, is shown for a Geminga-like pulsar 
with a dipole moment inclined at 22.5° to the rotation axis. The 
pulsar becomes significantly less observable at higher frequencies. 
Similar results are obtained for pulsars with different model pa- 
rameters. The axes show the base-10 logarithm of the photon 
energy in MeV along the horizontal, and the value of f2 p along 
the vertical. 



Setting £ = 36 (orthogonal fields) and writing n in 
terms of the Goldreich-Julian number density tigj, the 
gamma-ray efficiency is 



83- 



V-y = 



-Prot + Pi 



1+83- 



(27) 



Here, the number density has been written as a frac- 
tion of the Goldreich-Julian number density; the gamma- 
ray efficiency only depends on this ratio, not on the pulsar 
properties directly. 

As n/naj approaches unity, so does the gamma-ray ef- 
ficiency. This is seen in the cases of Geminga and PSR1055- 
52, where the observations suggest that ~ 1 (Mayer- 
Hasselwander et al 1994, Fierro et al 1993). 



6.3 The Distribution in Viewing Angle 

The emission is not beamed equally into all directions. Fig- 
ure ^ shows the emission (averaged over a period) at 1 MeV 
as a function of the viewing angle (9 V — to 90°) for two 
Vela-like pulsars: one with orthogonal fields (solid line) and 
the other with a dipole inclination of 11.25° (dashed line). 
The power is scaled to unity at the maximum. 

The fraction of the total solid angle into which the 



power is beamed can be quantified in terms of a power- 
weighted average: 

P(Q) 



(28) 



where Q represents solid angle, not the angular velocity of 
the pulsar. 

Figure ^ shows Q p as a function of frequency for a 
Geminga-like pulsar with a dipole moment inclined at 22.5° 
to the rotation axis. A general trend to lower Q p (and there- 
fore less observable pulsars) is seen for increasing frequency, 
with a much faster dropoff above the break in the spectrum 
(near several hundred MeV). Similar results are obtained for 
pulsars with different model parameters. 



7 SPECTRA OF EMITTED RADIATION 

Since curvature radiation is the only contribution to the light 
curves, the frequency distribution of curvature emission by 
a single particle moving instantaneously in circular motion 
must be considered (see, for example, Jackson (1975)). The 
shape of the curve depends on only two parameters - the 
energy of maximum intensity (E m = 0.14?i7 3 c/p, where 7 
is the charge's Lorentz factor, c is the speed of light, and 
p is the radius of curvature of the particle's path), and the 
maximum intensity. For energies much less than Em, the 
number distribution in frequency is proportional to _E~ 2 / 3 ; 
above Em, the curve drops off exponentially, and essentially 
disappears for E > 30i5 m . This function is plotted in Fig- 
ure m 

At low frequencies, therefore, each part of the light 
curve will fall off in photon number like E' 2 ^ 3 , and the 
shape of the light curve will remain essentially constant. At 
higher frequencies, the contributions from particles radiat- 
ing with different i5 m 's will result in a light curve which is 
variable with frequency. If several different locations in the 
magnetosphere make significant contributions to the power, 
the shape of the spectrum can deviate significantly from the 
single-particle spectrum. 

Figure ^ shows several different time-averaged spectra 
from model pulsars, all with /i — 10 31 G cm 3 , x — 90°, and 
seen from a viewing angle of 83.4°. Spectra for pulsars with 
three different angular velocities are shown: Q — 25 rad/s, 
75 rad/s, and 150 rad/s. The slow pulsar shows a very hard 
spectrum with an energy break near several GeV, but rel- 
atively featureless otherwise. The 75 rad/s pulsar shows a 
slight softening of the spectrum near a few hundred MeV, 
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Figure 8. Representative spectra for three model pulsars, all with fj, = 10 31 G cm 3 , \ = 90°, and seen from a viewing angle of 84.3°. 
The solid line corresponds to a pulsar with f2 = 150 rad/s, the dashed-dotted line to one with Q = 75 rad/s, and the dashed line to 
H = 25 rad/s. The faster pulsars show more spectral detail at low frequency. 



Figure 7. The curvature spectrum for a single particle in instan- 
taneously circular motion. Logarithm of photon energy in units 
of the reference energy h^c/p is shown along the horizontal, and 
logarithm of arbitrary differential photon number is plotted along 
the vertical. 



and then returning to a hard profile until an energy break 
near 100 GeV. The 150 rad/s pulsar has a significantly 
harder spectrum with no break even through the TeV re- 
gion. It is interesting to note that the photon flux for all 
three pulsars is relatively constant in the 100 MeV range 
even though the total power outputs vary considerably. 

Figure ^| shows time-averaged spectra for an orthogo- 
nal field, Geminga-like pulsar with a dipole moment of 10 30 
G cm 3 at a viewing direction of V — 84.4°. The spectrum 
does change slightly with viewing angle, though the ma- 



Figure 9. The model spectrum of a Geminga-like pulsar from 1 
MeV to 10 GeV. The data shown are EGRET data in the range 
40 MeV to 6.3 GeV, from Mayer-Hasselwander et al 1994. The 
theoretical curve has been shifted vertically on the plot to best fit 
the data. The horizontal axis shows the logarithm of the energy 
in MeV, and the vertical axis shows the logarithm of differential 
photon number in units of photons/s/cm 2 /MeV. 



jor features (location of the energy break and the shape 
of the spectrum) are roughly the same. Data from Mayer- 
Hasselwander et al 1994 are shown in the interval 40 MeV 
to 6.3 GeV; the theoretical curve has been shifted vertically 
on the plot to best fit the data. 

The different regions of the light curve can have some- 
what different spectral properties. Figure [To| shows the light 
curve (c) from Figure ^ separated into four components. The 
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Figure 10. The frequency dependence of the light curve (c) from Figure 4. (a) shows the light curve at 1 MeV, and (b) shows the 
emission at 10 GeV. (c) shows the spectra for four different components in the light curve. The line dashing of the components in the 
light curves corresponds to that of the lines in the spectral plot. The axes show the base-10 logarithm of energy in MeV along the 
horizontal and logarithm of arbitrary differential photon number along the vertical. 



components correspond roughly to the peaks at low energy 
(1 MeV), shown in (a), and the high-energy (10 GeV) peaks, 
which fall at phases of zero and 0.5. 

The phase of the gamma-ray peak tends to change with 
energy when the energies are close to the energy break if 
it is significantly different from the phase of the radio peak. 
This occurs because the emission of highest frequency comes 
from particles moving nearest to the star, where the fields 
are larger. As the energy increases past the break, the peaks 
tend to move into phase with the radio peaks, since both are 
coming from relatively close to the star near the poles. This 
alignment of gamma-ray phase with radio phase at high fre- 
quencies is a provable prediction of the model, given the pre- 
sumed location of radio emission. It is difficult to see this ef- 
fect in the present observations, since the phase change only 
happens at very high frequencies (typically above 10 GeV). 
However, if the low-energy gamma peaks are not aligned 
with the radio peaks, it does lead to the "interpulse" re- 
gions at low energy having harder spectra and higher energy 
breaks than the low-energy pulses, as is seen, for example, 
in the case of the Vela pulsar (Kanbach et al 1994). 

Generally the low-energy spectral index of —2/3 is 
harder than observed in the X-ray region for those pulsars 



where this has been measured. However, this is precisely the 
region where the model begins to break down, and radiation 
from other areas (such as the inner magnetosphere, where 
unperturbed curvature radiation likely is not the dominant 
contribution to the spectrum) becomes important. 



8 DISCUSSION 

A new model of high-energy pulsar emission has been pre- 
sented, with the simplifying assumption that the charge den- 
sities and current densities are small enough that the elec- 
tromagnetic fields can be represented by the Deutsch fields 
in the outer magnetosphere. 

The assumption that the magnetosphere can remain 
evacuated is justified by recognising that charges created 
near the star follow paths bound to the inner magnetosphere, 
and cannot short out the fields farther out. We create an ap- 
proximate self-consistent global solution to the electromag- 
netic fields and charge motions by extending the "stellar" 
radius (inside which the fields are frozen-in) to 0.1 of the 
light cylinder radius; charge created in the inner magneto- 
sphere can flow out along null surfaces where the electric 
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field is perpendicular to the magnetic field. The outer mag- 
netosphere is the location of particle acceleration and high- 
energy emission. 

An important approximation has been developed to 
chart particle motion in fields where the timescale required 
to adjust to the field direction is much less than the timescale 
to radiate away the charge's momentum perpendicular to 
the magnetic field, called the drift frame bead-on-a-wire 
(DFB) approximation. This reduces the problem of parti- 
cle motion to a simple first order, ordinary set of differen- 
tial equations, reducing computation time enormously. This 
approximation could prove useful for developing global, self- 
consistent models of the pulsar magnetosphere, as it gives a 
relatively simple description of charge paths in terms of the 
electric and magnetic fields. 

Several constraints on the applicability of the model can 
be made, the most important of which is that the number 
density required to produce the observed power output must 
be much less than the Goldreich- Julian charge density. All 
observed gamma-ray pulsars satisfy this constraint, as well 
as that which maintains the dominance of curvature radi- 
ation to the high-energy spectrum in the model framework 
and that which maintains the applicability of the DFB ap- 
proximation. 

The light curves generated from the model match the 
general properties of the observed light curves quite well. 
When the magnetic dipole moment is orthogonal to the rota- 
tion axis, the light curves are mainly double-peaked; the sec- 
ond half of the light curve is identical to the first. Other incli- 
nations can give single peaks, non-symmetric double peaks, 
and triple peaks, as well as practically unpulsed emission in 
some cases. 

The biggest weakness of the model is its present in- 
ability to predict the actual power output of the pulsars in 
gamma-rays: it does not assume any value for the number 
density at the starting sphere, only that it be much less 
than the Goldreich- Julian number density so that the elec- 
tric fields are not shorted out. It also says nothing about 
any radiation produced in the inner, bound magnetosphere, 
where particle energies will be somewhat lower than those 
in the outer, evacuated magnetosphere. This might provide 
significant emission at lower energies. 

Our model compares quite favourably with the other 
popular models for pulsar gamma-ray emission. The spec- 
tra predicted by our model match the observations in some 
cases better than those generated by the outer gap model 
of Chiang & Romani (1994), and as well as those presented 
for the extended polar cap model of Daugherty & Harding 
(1995). Our model also has the potential to generate more 
complex light curves than either the outer gap or polar cap 
models, which can only produce single or double peaks. How- 
ever, this is something of a double-edged sword, as a random 
sample of light curves would be somewhat less regular than 
is observed. 

This work has introduced the general Deutsch field 
gamma-ray pulsar model; applications of the model to the 
known gamma-ray pulsars are developed in Paper II, which 
also considers several statistical properties of a population 
of these pulsars. 
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